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INTRODUCTION 


Numerous models of physical systems contain parameters whose values are 
not known exactly. This paper attempts to address some of the physical and 
mathematical complexities arrising in the prediction of. the statistical be- 
havior of such systems. Although the discussions in the paper are far from 
providing a satisfactory solution to such problems, they perhaps, by utiliza- 
tion of simple examples, will create a greater awareness of the statistical 
effect of random parameters. 


PROBLEM FORMULATION 


Consider the problem of determining the statistical properties of the re- 
sponse of a finite dimensional linear dynamical system with random coefficients 
(constant with respect to time) and subjected to stochastic forces. Mathe- 
matically the problem is represented by the following equation: 


= Ax(t) + f(t) 0<t 

x(o) = X 

o 


( 1 ) 


where x(t), f(t), and x are n~dimensional random vectors, A is an nxn random 
matrix. The problem is to determine statistical properties (mean value, 
variance, correlation fimction, spectral density, distribution, etc.) of x(t) 
knowing the statistical properties of x^, f(t), and A. 

EXISTENCE OF SOLUTION 


If the derivative in equation (1) is interpreted in the almost sure sense 
then existence and uniqueness of a solution follows from appropriate results 
in and if f(t) is almost surely continuous then a solution in this sense 
would exist and be given by: 
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f (x)dT 


( 2 ) 


x(t) = 

o o 


gACt-x) 


Howeve^ since the discussions in this paper will be concerned with second 
moments of the solution, it would seem appropriate to require the derivative 
in (1) to be a mean square derivative and to consider the differential equa- 
tion (1) over the Hilbert space where Z denotes the space of second order 
random variables. If A is a bounded operator over Z*^ (probably equivalent to 
requiring A to be almost surely bounded) then the theory of ordinary differ- 
ential equations would yield a unique solution given by (2), where the inte- 
gral was a mean square integral, provided f(t) is mean square continuous. In 
general, however, A may not be bounded, l.e. the product of two second order 
random variables will not be a second order random variable and then the 
appropriate theory discussing existence of a solution to (1) would likely be 
a requirement that the solution be the action of a semigroup on the initial 
condition. For example, if there exists a real number such that: 

I 1a[A+A^ “ a] ^11 ^ C, for all complex A with Re A>0, a real 


number and [ | | | denoting the norm over r’^, then there will be a solution in 
the mean square sense to (1) and further: 

E [ (e^^x^ ) *^(e^^x^ ) ] <_ C2e ° E [x^ ^x^ ] 

In particular if A can be chosen to be negative then the solution will be 
asymptotically stable. This approach to the problem exhibits a solution with 
the only stipulations that x eZ’^ and f (t) be mean square continuous. Another 
approach to finding a mean square solution to (1) would be to require condi- 
tions on Xq, a, f(t) such that (2) is a solution to (1). If x , A, and f(t) 
are mutually independent, then requiring e^^ and Ae^^ to have second moments 
would insure that (2) satisfies (1) . The following elementary examples attempt 
to illustrate the above discussion. 


EXAl^LES 


Example 1 

Consider the first order homogeneous equation (n = 1) . 


dx 

dt 


ax 


x(o) = X 
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with a uniformly distributed between a and 3. Clearly a is a bounded operator 
over Z thus, for example, if x^ is independent of a it follows that: 


'■’‘o' 

If 3>0 then E[x(t)] becomes arbitrarily large even if the mean of a is nega- 
tive. 

Example 2 

Consider the 


P^(a> - 


Although a is not 

for all Re X>0. Thus a solution exists. If x^ is independent of a then it 
can be sho^ that 

E[K(t)]=g^ Elx^]. 

Note again that if 3>0 E[x(t)] becomes arbitrarily large. 

! 

Example 3 

Consider the above example with a Gaussian with mean u and variance o. 
Clearly a is not bounded and further no Xq can be chosen to make i X 

X+X —a 

at at ^ 

However if a is independent of x then ae and e do have second moments and 

it follows that; ° 


1 C. 


P 

above problem with the density of a given by; 4 



a bounded operator over Z clearly tT; < C if X >3 

X+X -a — o 


E[x(t)] = exp } E[x ] . 

I o 

a 

However, regardless of a and p, E[x(t)3 becomes arbitrarily large. 
Example 4 

Consider the above example with the density of a given by: 
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Y 






strated that x(t) does not have a first moment for t>a. Thus it makes no 
sense in this problem to attempt to calculate E[x(t)]. 


STATIONARY RESPONSE AND SPECTRAL DENSITY 


Assume that the existence of the solution in the mean square sense to (1) 
is known and is expressible as: 

/^\ At ^ -t A(t-T) j;/ SJ 
x(t) = e X t / e f(T)dT 

o o 

If A, X , and f(t) are mutually independent and further if f (t) is stationary 
with correlation matrix R^ then it follows that: 

E[(x(t+A) - E[x(t+A)]) (x(t) - E[x(t)])^] = 

T 

Rf ^ ] dn^ <1^2 

If it can further be shown that ||e^^|| Ce^^ with 8<0, then it follows in 
the usual way that as t goes to “ x(t) becomes stationary with: 

By taking the Fourier transform of Rj(A) it is easily shown that the spectral 
density of x(t) is given by: 

S (w) = E[ lA+io)] ^ S^Cm) [A-lai] ^ ] 

X r 
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Example 5 


Of the previous examples only Example 1 with 3<0 and Example 2 with 3<0 
eventually have stationary solutions . In example 1 (with 3<0) it is easily 
seen that : 


S (o)) = ^ [tan ^ tan ^ ~ ] S^(o)) 
X 0) CO 0) •' f 



If a was not a random variable then S (to) == ” 9 ^ "y 

a +(o 


and a plot of this follows: 



Those readers who have gotten to this point in the paper recognize 
it as a fraud. The paper (1) presents a physical problem, i.e.: how do you 
calculate the statistical properties of the response of dynamical systems 
which have random parameters, ( 2 ) presents possible mathematical models that 
pertain to the physical problem and (3) presents, via simple examples, where 
the problems are in trying to solve the problem. The result in example 3, 
where a is Gaussian, shows that regardless of how negative the mean value and 
how small the variance of a, the mean value of the solution goes to time 

goes to <». In particular, it makes no sense to talk about the spectral density 
of the solution. 

In the opinion of the author closed form solutions to problems beyond n=l 
are not feasible and current work centers around the study of the accuracy of 
approximate methods that have been proposed in the literature. 
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